ADA108844 


SAMSO  TR-79-9 


SAI  DOCUMENT  NO.  SAI-79-505-VF 

PERFORMANCE  TECHNOLOGY  PROGRAM 
j  (PTP-S  II) 

VOLUME  II 

A  THREE-DIMENSIONAL  INVISCID  BLUNT  BODY  CODE 
FOR  ASYMMETRIC  NOSETIPS  IN  PITCH  AND  YAW 

SCIENCE  APPLICATIONS,  INC. 

APPLIED  MECHANICS  OPERATION 
WAYNE,  PENNSYLVANIA  19087 

DECEMBER,  1978 


FINAL  REPORT  FOR  PERIOD  MARCH  1978  -  DECEMBER  1978 


CONTRACT  NO.  F04701-77-C-01 26 


APPROVED  FOR  PUBLIC  RELEASE; 


AIR  FORCE  SPACE  AND  MISSILE 
SYSTEMS  ORGANIZATION 
LOS  ANGELES,  CALIFORNIA 


DISTRIBUTION  UNLIMITED. 


DTIC 

ELECTS 
DEC  38  1881 


81  12  23  il8 


’a?'JT?w  •  7*^*1 

V"" 


This  final  report  was  submitted  by  Science  Applications, 
Inc.,  1200  Prospect  Street,  La  Jolla,  California  92037,  under 
contract  number  F04701-77-C-0126  with  the  Space  and  Missile 
Systems  Organization,  AFSC,  Los  Angeles,  California  90009* 
Capt.  Robert  J.  Chambers,  SAMSO/DYSE  was  the  Project  Officer 
in  charge.  This  technical  report  has  been  reviewed  and  is 
approved  for  publication. 


— • 


CORMACK,  LT  COL,  DSAP 
Reentry  Technology  Division 


FOR  THE  COMMANDER 


WlLLIAM-GQL»BfeG,  COL,’  \&fjT 
Director,  Systems  Technology 


A 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PACE  (When  Do la  Entered) 


REPORT  DOCUMENTATION  PAGE 


I.  REPORT  NUMBER 


KEAD  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


i.  recipient's  catalog  number 


SAMSO  TR-  79-9 


4.  TITLE  (end  Submit  I 

Performance  Technology  Program  (PTP-S  II) , 
Vol.  II  -  A  Three-Dimensional  Inviscid  Blunt 
Body  Code  for  Asymmetric  Nosetips  in  Pitch 
and  Yaw 


7.  AUTHOR/*; 

Darryl  W.  Hall 
Charles  L.  Kyriss 
Catherine  M.  Dougherty 


S.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Science  Applications,  Inc. 

994  Old  Eagle  School  Rd . ,  Ste.  1018 
Wayne,  Pennsylvania  19087 


II.  controlling  office  name  ano  AODRF«« 

Space  and  Missile  Systems  Organization 
(SAMSO) 

Los  Angeles,  California 


4.  MONITORING  AGENCY  name  4  AODRE5S/I/  dlllerenl  Iron i  Controlling  Olllcm) 


w 


S.  TYPE  OF  REPORT  t  PERIOD  COVERED 

Final  Report:  3/78-12/78 


S.  PERFORMING  ORO.  REPORT  NUMBER 

SAI-79-505-VF 


t.  contract  or  grant  NUMBER/*) 

F04701-77-C-0126 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  •  WORK  UNIT  NUMBERS 


Task  3. 2. 1.1 


12.  REPORT  date 

December  1978 


IS.  NUMBER  OF  PAGES 

22 


IS.  SECURITY  CLASS,  (ol  Ihlo  report) 


UNCLASSIFIED 


Sa.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


16-  DISTRIBUTION  STATEMENT  (ol  thlt  Ropott) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (o  t  Iko  obolroei  on  to  rod  In  Block  20,  II  dllloront  tram  Report) 


IS.  KEY  WORDS  (Continue  on  rovoroo  o Ido  II  noceeoery  end  Identity  by  bloc k  number) 

Flow  Field 
Blunt  Body  Code 
Reentry  Vehicles 
Asymmetric  Nosetips 

/  _  _ _ 


10/  ABSTRACT  (Cottllnuo  on  r«v«r««  aid*  If  nac mmoory  *nd  Identify  by  block  numbrnr) 

/ 

/.This  report  documents  modifications  made  to  the  Moretti  inviscid 
three-dimensional  blunt  body  code  to  remove  the  assumption  of  a 
pitch  plane  of  symmetry.  Details  are  provided  on  the  analysis 
required  to  extend  this  code  to  treat  fully  asymmetric  nosetip 
geometries  at  both  angle  of  attack  (a)^and  sideslip  10/  Modi¬ 
fications  to  the  code  inputs  required  for  this  capability  are  also 
described.  ^ 


DD  1473 

IF  A  CBIMI  LEI 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Wh*n  Dele  Entered) 


n  onknijk  Lo+rnl  A  AWfa.  A  .1  ojo  ,  o.~... 


TABLE  OP  CONTENTS 


i 


ABSTRACT  (DD1473)  .... 
TABLE  OF  CONTENTS  .... 
SECTION  1  INTRODUCTION.  .... 

SECTION  2  ANALYSIS  ..... 

2.1  Governing  Equations  at  Centerline 

2.2  Treatment  of  Sideslip 

2.3  Derivative  Approximations 

2.4  Body  Geometry 

2.5  Modifications  to  Code  Inputs 

SECTION  3  VALIDATION  ..... 

SECTION  4  CONCLUSIONS  . 

SECTION  5  REFERENCES  ..... 


1 

2 

4 

4 

10 

11 

12 

13 

14 
16 
17 


SECTION  1 


INTRODUCTION 


This  report  documents  modifications  made  to  the  three- 
dimensional  inviscid  Moretti  blunt  body  flow  field  code  to 
remove  the  assumption  of  a  pitch  plane  of  symmetry.  These 
modifications  make  possible  calculations  at  sideslip  (6)  as 
well  as  angle  of  attack  (a) ,  and  also  allow  calculations  to  be 
made  on  asymmetric  nosetip  geometries  that  do  not  have  a  pitch 
plane  of  symmetry. 

The  blunt  body  code  being  modified  is  based  on  the 

analysis  of  Moretti  and  was  first  documented  by  Kyriss  and 

(2) 

Harris  .  A  more  detailed  exposition  of  this  original  analysis 
is  given  in  Reference  3,  The  ability  of  the  Moretti  blunt  body 
code  to  accurately  compute  inviscid  flows  over  ablated  nosetip 
geometries  (subject  to  the  restrictions  of  the  assumption  of 
a  pitch  plane  of  symmetry)  has  been  investigated  and  demon- 

14) 

strated  by  Hall,  Kyriss,  Truncellito  and  Martellucci  ' . 

The  Moretti  blunt  body  code  is  based  upon  an  inviscid, 
time-dependent  technique  which  solves  the  governing  equations 
in  non-conservation  form  using  the  second-order  accurate 
MacCormack ^  finite  difference  scheme.  The  equations  are 
formulated  in  a  spherical-polar  coordinate  system,  and  either 
ideal  or  equilibrium  air  thermodynamics  may  be  used  in  the 
calculations . 

Special  forms  of  the  governing  equations  are  required 
along  the  singular  centerline  of  the  spherical-polar  coordinate 
system.  In  the  original  form  of  this  blunt  body  code,  these 
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special  forms  were  derived  under  the  assumption  of  a  pitch 
plane  of  symmetry,  in  which  many  of  the  terms  in  the  pitch 
plane  vanish. 

As  part  of  the  current  effort  to  remove  the  pitch  plane 
of  symmetry  assumption,  special  limiting  forms  of  the  governing 
equations  have  been  developed  to  allow  treatment  of  non- 
symmetric  geometries  and  sideslip.  The  new  centerline  equations 
are  presented  in  Section  2.1  of  this  report.  The  treatment  of 
sideslip  at  shock  points  is  discussed  in  Section  2.2,  and  the 
finite  difference  expressions  required  by  the  new  centerline 
equations  are  described  in  Section  2.3.  Changes  made  to  the 
body  geometry  specification  procedure  are  given  in  Section  2.4. 
Finally,  changes  in  the  code  inputs  required  by  this  modifi¬ 
cation  effort  are  given  in  Section  2.5. 


2.1  Governing  Equations  at  Centerline 

In  the  (r,8,$)  spherical-polar  coordinate  system, 
shown  in  Figure  1,  the  governing  inviscid  equations  may  be 
written  as: 

Pt  +  uPr  +  vPQ/r  +  wP<j,/r  sin  9  +  V(ur  +  vQ/r 

+  w,/r  sin  8  +  2u/r  +  v  cot  Q/r)  *  0  (1) 

f 

ufc  +  uur  +  vUg/r  +  wu^/r  sin  9  -  (v2  +  w2)/r 

+  pP„/p  *  0  (2) 


vfc  +  uvr  +  Wg/r  +  wv^/r  sin  9  +  uv/r 

-  w2  cot  0/r  +  pPg/pr  =  0  (3) 

I 

7 

i 

wfc  +  uwr  +  vwg/r  +  ww^/r  sin  9  +  uw/r 

4-  vw  cot  8/r  +  pP^/pr  sin  9=0  (4)  \ 

1 

i 

s  +  us,  +  vsfl/r  +  ws./r  sin  9=0,  (5)  j 

c.  r  a  p  ; 

I 

I 

where  P  is  the  logarithm  of  pressure,  u,v,  and  w  are  the  r,0, 

and  <j>  velocity  components,  respectively,  and  3  is  the  entropy.  ■ 


I 

I 

ii*h  f  | 
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At  the  centerline,  where  9  -  it,  sin  8  »  0  and  cot  0  *  * 
these  equations  are  singular.  Removing  the  singularities 
results  in  the  following  special  forms  of  the  equations,  valid 


only  along  the  centerline: 

Pt  +  uPr  +  vPg/r  -  wP^/r  +  y(ur  +  2vQ/r 

+  2u/r  -  wQ<{,/r)  *  0  (6) 

ufc  +  uur  +  vUg/r  -  wug^/r  -  (v2  +  w2)/r 

+  pPr/p  =0  (7) 

vfc  +  uvr  +  vvQ/r  -  wvg^/r  -  wwQ/r  +  uv/r 

+  pPg/pr  *  0  (0) 

wt  +  uwr  +  vw0/r  "  we^r  +  w0/r 

+  uw/r  -  pPg^/pr  *  0  19) 

st  +  usr  +  vsg/r  -  wsg^/r  =  0  .  (10) 


In  the  previous  analysis,  with  the  assumption  of  a 
pitch  plane  of  symmetry,  these  equations  are  used  only  in  the 
plane  of  symmetry,  where  w  =  0,  wQ  =  wr  =  0,  and  Pg^  *  Sg^  * 
Ug^  =  vg^  =  0,  thus;  simplifying  the  analysis. 

The  generalized  equations  valid  at  the  centerline. 
Equations  (6)  -  (10) ,  are  transformed  to  the  computational 
space  using  the  transformation  defined  by: 

Z  =  [r-rb(9,4>)  ]/[rs(0,4>,t)  -  rfa (0 , <J>)  ] 

Y  =  TT  -  9 

X  =  4> 


(11) 

(12) 

(13) 
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where  the  body  surface  is  defined  by  r  *  r^(8,$)  and  the  bow 
shock  surface  by  r  ■  ra(9,$,t).  Derivatives  are  transformed 

0 

according  to  the  chain  rule,  resulting  in 

a  1  a 


3r  “  Z  Zz 

(14) 

T5"  Ty  +  Z  t(Z-l) rb^-ZrSQ] 

a 

Zz 

(15) 

=  lx  +  J[(z"1)rb^  “Zrs$l 

a 

Zz 

(16) 

3  _  3  Zicf  3 

Tt  ~  Zr  ~  — —  Zz 

(17) 

where  5  =  r  (0,4>,t)  -  rb(9,<J>). 

Also,  at  the  centerline: 

aoa?  =  ~  txaY  +  I  t  iz 

+  |  t(Z-l)rbe-Zrs&]  .  (18) 

Following  the  coordinate  transformation,  the  governing 
equations,  both  on  and  off  the  centerline,  may  be  written  in 
one  form,  where  certain  coefficients  take  different  values  at 
the  centerline  from  their  values  elsewhere.  The  resulting 
general  forms  of  the  equations  may  be  written  as 

PT  +  APy  +  NPX  +  FPZ  +  Mx  +  Y[uz/5  +  Bvz 

+  Ewx  +  Jwz  +  W  +  (2u-vy)/r  +  H]  =  0  (19) 


6. 


(20) 


uT  +  Auy  +  Nux  +  FUj  +  M2  -  (v*  +  W*)/r 
+  pP2/«P  -  0 

vT  +  Avy  +  Nvx  +  Fvz  +  M3  -  Au  -  K 
4“  p(BPz  -  PY/r)/P  *  0 


wT  4-  Awy  +  NWX  ♦  Fwa  +  M4  +  uw/r  ♦  L 
+  p(EPx  +  JP2  *  F)/P  "  0 

sT  +  A3y  +  NSX  ♦  FSZ  +  M5  -  o 


(22) 

(23) 


Coefficients  whose  expressions  ere  the  same  both  on  and 
off  the  centerline  are 


C  *  (Z-l)rbQ-ZrSQ 

q  *  (z-1) rb^-2rs^ 

A  »  -v/r 
B  s*  C/rd 
N  <  Ew 

t  *  (u-2rSt  -  AC  4-  NQ)/6 
J  *  EQ/d 


jT. T  _  ' 


At  points  off  the  centerline,  the  remaining  coefficients  take 
the  forms 


B  ■  1/r  sin  9 
H  *  Ev  cos  0 
K  ■  Nw  cos  9 
L  ■  v»H 
P  »  0 
W  »  0 

M1  *  M2  =  M3  =  M4  50  M5  “  0 

On  the  centerline  these  coefficients  are  expressed  as 
E  -  0 

H  =  (-vy  +  Cvz/«)/r 
K  =  2w(-wy  +  Cwz/6)/r 
L  =  wH  +  v  (-wy  +  Cwz/6)/r 

P  =  -(-PXY  +  CpXZ/6  +  pzQ9/6)/r 
W  =  - (-wXY  +  Cwxz/5  +  wzQQ/6)/r 
=  wP 

Mj  =  -w(-Uxy  +  +  UZ®0^^^r 

M3  =  -w  ( -vxy  +  Cvxz/S  +  vZQ0/6)/r  +  w^~wy  +  Cwz/6)/r 

M4  =  wW  -v(-wy  +  Cwz/6)/r 

M5  =  -w(-sxy  +  Csxz/5  +  szQQ/6)/r 


8 


where 


Qe  *  ze  (rtxj>”rs4>>  +  (z-i)rb0(j>  “  Zra 

Along  the  centerline,  calculations  are  performed  only 
in  the  <$>  *  0  plane.  The  flow  variables  P,s,  and  u  at  the 
centerline  are  the  same  for  all  values  of  <fc ;  the  v  and  w 
velocity  components  for  other  values  of  $  may  be  expressed  as 


v(<J>)  «  v (0) 

cos  <}> 

-  w(0) 

sin  $ 

(24) 

w(4>)  *  v(0) 

sin  <{> 

+  w(0) 

cos  $ 

• 

(25) 

The  equations  presented  above  are  used  directly  at 
field  points  on  the  centerline.  At  body  and  shock  points,  the 
computational  procedures  used  (which  are  described  in  References 
2  and  3)  require  a  linear  combination  of  these  equations;  thus 
the  body  and  shock  point  procedures  on  the  centerline  require 
modification  to  eliminate  the  plane  of  symmetry  assumption. 

These  modifications  follow  directly  from  Equations  (19)  -  (23)  . 

Without  a  pitch  plane  of  symmetry,  however,  the  ex¬ 
pressions  for  the  body  and  shock  normals  at  the  centerline 
must  be  modified  to  allow  for  components  normal  to  the  pitch 
plane,  which  previously  were  zero.  The  outward  body  normal 
at  the  centerline  is  now  expressed  as 

er  -  rb9/rb  e9  +  rb9lf/rb  8^ 

nb  "  |1  +~b0*/rb»  +  rb04,i/rb*!  V*  (26) 

and  the  inward  shock  normal  becomes 


9. 


(27) 


A 


~er  :  rae/r3  e9 
[1  +  rS9Vr3!  + 


-  rs9</r3  •* 

rs8*!/r^,1/! 


2.2  Treatment  of  Sideslip 


The  addition  of  the  capability  to  handle  sideslip 
(3  7*  0)  requires  only  that  the  definition  of  the  freestream 
velocity  vector  be  altered.  Defining  the  angle  of  attack  (a) 
and  the  sideslip  angle  (3)  as  shown  in  Figure  2,  the  freestream 
velocity  vector  may  be  written  in  terms  of  the  (r,0,<{>)  spheri¬ 
cal-polar  coordinate  system  as 


ue  +  vm  ea 
oo  r  co  y 


+  w 

00 


A 


where 


uoo  =  qootGOS  a  cos  9  cos  9  +  (sin  B  sin  <j> 

+  sin  a  cos  &  cos  $  ) sin0) 

v  =  [-cos  a  cos  3  sin  0  +  (sin  3  sin  <j> 

OO  CO 

+  sin  a  cos  6  cos  <J>  )  cos  0  ] 


w  =  q  [sin  3  cos  <b  -  sin  a  cos  3  sin  d>  ] 

OO  *00 


with  qw  =  |vj 


(28) 


(29) 


(30) 

(31) 


The  presence  of  a  non-zero  3  also  effects  the  definition 
of  the  equivalent  body  angle  used  to  initialize  the  surface 
pressures  and  shock  shape.  The  equivalent  body  angle  is  now 
defined  as 


0bEQ 


0fa  -  a  cos  $  -  3  sin  <J> 
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(32) 


where 


(33) 


-(rbg/r^)  sin  0  -  cos  6 
-(rbg/rk)  cos  6  +  sin  0 

Approximations 

The  generalization  of  the  limiting  forms  of  the  govern¬ 
ing  equations  at  the  centerline  introduces  second  derivatives 
that  previously  vanished  in  the  pitch  plane  of  symmetry.  The 
finite  difference  approximations  to  these  derivatives  may  be 
written  as 

s 2  f  —  —  _  _ 

=  [f  (N,M  +  1,  L  +  1)  -  f  (N/M,  L  +  1) 

-f (N,M  +  1 / L)  +  f (N,M,L) ]/(AX) (AY)  (34) 

and 

-  [f  (N  +  1,  M,L  +  1)  -  f  (Nr  Mr  L  +  1) 

-f(N  +  1,  M,L)  +  f  (N,M,L)  ]/(AX)  (AZ)  (35) 


0^  =  tan-1 

2.3  Derivative 


where  Z  =  N  AZ 
Y  =  M  AY 
X  =  L  AX 

At  the  centerline  the  barred  indices  take  the  forms 

N  =  maxflr  min{N-I  +  lr  NMAX  -  1})  (36) 

M  =  1  (37) 

L  =  1  (I  =  1) 

LMAX  -  1  (1=2)  (38) 


11 


where  I  =  1  in  the  predictor  stage  and  I  =  2  in  the  corrector 
stage  of  the  MacCormack  finite  difference  scheme.  NMAX  and 
LMAX  are  maximum  values  of  the  Z  and  X  indices,  respectively. 
Note  that  without  a  pitch  plane  of  symmetry,  L  =  1  (<)>  =  0°) 
and  L  =  LMAX  (<J>  ®  360°)  correspond  to  the  same  physical  plane 
(periodicity  condition) . 

The  calculation  of  r^g^  and  rgg^  is  performed  using 
Equation  (34),  noting  that  (from  Equation  (18)) 

g  2  g  2 

=  “  9X3Y  *  (39) 

2.4  Body  Geometry 

In  the  original  blunt  body  code,  two  options  exist  for 
the  definition  of  the  body  geometry,  as  described  in  Reference 
6.  The  first  geometry  option  requires  the  specification  of 
geometric  profiles  in  a  number  of  <J>  planes;  these  profiles  may 
be  defined  either  analytically  or  by  tabular  input.  The  second 
option,  for  "bi-elliptic"  geometries,  requires  the  definition 
of  profiles  in  the  "lee"  (<j>  =  0°)  and  "wind"  (<j>  =  180°)  planes, 
and  optionally,  the  "side"  plane,  by  tabular  input. 

The  extension  of  this  code  to  remove  the  pitch  plane 
of  symmetry  includes  the  capability  of  defining  geometric 
profiles  (either  analytically  or  by  tabular  input)  for  a  range 
of  <})  from  0°  to  360°.  No  modifications  to  the  input  variables 
have  been  made.  However,  if  a  fully  asymmetric  geometry  is  to 
be  specified,  profiles  must  be  defined  in  both  the  =  0°  and 
<}>  =  360°  planes  (which  are  identical)  . 
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SECTION  3 


VALIDATION 


The  validation  of  the  blunt  body  code  modifications 
developed  in  this  effort  to  remove  the  pitch  plane  of  symmetry 
assumption  has  been  carried  out  by  comparisons  of  numerical 
results  obtained  with  both  the  original  and  modified  codes. 

This  validation  procedure  is  possible  because  the  modifica¬ 
tions  required  to  remove  the  assumption  of  a  pitch  plane  of 
symmetry  have  not  altered  the  fundamental  numerical  algorithm, 
and  the  original  blunt  body  code  has  been  previously  vali¬ 
dated  through  extensive  comparisons  of  numerical  results  to 
experimental  data,  as,  for  example,  in  Reference  4. 

Because  the  basic  numerical  algorithm  has  not  been 
altered,  the  modified  blunt  body  code  developed  in  this  effort 
is  subject  to  the  same  limitations  on  its  applicability  as 
the  original  code,  except  that  flows  without  a  pitch  plane 
of  symmetry  can  now  be  treated.  These  limitations  have  been 
investigated  and  discussed  in  Reference  4. 

The  first  case  examined  in  the  validation  of  the 
modified  blunt  body  code  was  the  simple  case  of  a  sphere. 
Calculations  were  made  with  the  original  code  at  a  =  5°  {with 
a  pitch  plane  of  symmetry)  and  with  the  modified  code  at 
a  =  0°,  8=5°.  Comparisons  of  the  results,  rotated  to  account 
for  the  change  in  the  wind  vector,  revealed  no  significant 
differences  in  the  two  calculations.  As  a  further  step  in  the 
validation  process,  afterbody  solutions  were  computed  using 
initial  data  from  the  two  blunt  body  solutions,  and,  after 
the  appropriate  rotation  of  the  results,  no  significant  dif¬ 
ferences  were  found  in  the  afterbody  calculation  in  either 
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the  primary  flow  field  parameters  (e.g.,  pressure)  or  the 
integrated  forces  and  moments. 

Similar  validation  cases  were  run  with  nosetip  geom¬ 
etries  that  were  not  ax i symmetric,  but' which  did  have  pitch 
planes  of  symmetry.  (Calculations  at  sideslip  with  the  modi¬ 
fied  code  required  the  appropriate  rotation  of  the  nosetip 
geometry  so  that  the  plane  of  geometric  symmetry  would  be 
aligned  with  the  wind  vector.)  Again,  after  rotation  of  the 
results,  it  was  found  that  the  original  and  modified  codes 
produced  equivalent  results. 
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SECTION  4 


CONCLUSIONS 


An  existing  inviscid  blunt  body  code  has  been  modi¬ 
fied  to  remove  the  assumption  of  a  pitch  plane  of  symmetry, 
allowing  calculations  to  be  made  at  both  angle  of  attack  and 
sideslip  for  fully  asymmetric  nosetip  geometries.  All  capa¬ 
bilities  of  the  original  code  have  been  retained;  both  the 
original  and  modified  codes  are  still  subject  to  the  following 
limitations  on  their  applicability: 

1. )  The  shock  layer  flow  must  be  entirely 

inviscid,  except  for  a  thin  boundary 
layer  at  the  body  surface. 

2. )  The  nosetip  shape  must  permit  accurate 

flow  field  calculations  using  a  spheri¬ 
cal  coordinate  system. 

3 .  )  There  can  be  no  embedded  shocks  in  the 

transonic  flow  field. 

The  modifications  to  the  blunt  body  code  have  not 
increased  the  core  storage  requirements  above  those  of  the 
original  code,  and  there  has  been  no  change  in  the  computer 
time  requirements  (per  grid  point  per  iteration) . 

The  code  modifications  described  in  this  report  are 
available  in  an  UPDATE  format  for  the  SCOPE  or  N0S/BE1  operating 
systems.  Qualified  users  may  obtain  these  modifications 
through  SAMSO  (RSSE) . 
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